Sympathetic cooling of an atomic Bose-Fermi gas mixture 



W. Geist, L. You, and T.A.B. Kennedy 
School of Physics, Georgia Institute of Technology 
Atlanta, Georgia 30332-0430 
(February 1, 2008) 

Sympathetic cooling of an atomic Fermi gas by a Bose gas is studied by solution of the coupled 
quantum Boltzmann equations for the confined gas mixture. Results for equilibrium temperatures 
and relaxation dynamics are presented, and some simple models developed. Our study illustrate that 
a combination of sympathetic and forced evaporative cooling enables the Fermi gas to be cooled to 
the degenerate regime where quantum statistics, and mean field effects are important. The influence 
of mean field effects on the equilibrium spatial distributions is discussed qualitatively. 
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I. INTRODUCTION 



The magnetic trapping and forced evaporative cooling of alkali metal vapors has lead to the Bose-Einstein con- 
densation of a number of different atomic species [Q. As a consequence, the many-body physics of confined, weakly 
interacting Bose gases are now amenable to experimental investigation. To date most experimental and theoretical 
research on degenerate atomic gases has focused on bosons [pi . Recently, some exotic systems have been investigated 
experimentally, these include the Bose condensation of spin mixtures of 87 Rb using sympathetic cooling j^], and a 
report of spin domains of 23 Na in an optical trap Qj. There is now a growing interest in the properties of degenerate 
atomic Fermi gases and boson- fermion mixtures , although to date a degenerate atomic Fermi gas has not 

been achieved. 

In this paper we investigate the process of sympathetic cooling of an initially non-degenerate Fermi gas to the 
quantum degenerate regime using a Bose gas as coolant. Sympathetic cooling of a Fermi gas to degeneracy is necessary 
as a result of the suppression of s-wave scattering between identical fermions in spin symmetric states. Evaporative 
cooling of a pure fermion gas trapped in a single hyperfine state, which depends on rethermalization through atomic 
collisions, is thus ineffective at temperatures sufficiently low that only the lowest few partial waves contribute. The 
Fermi gas may instead be cooled by thermal contact with a cold Bose gas, which may be either already condensed 
before thermal contact with the Fermi gas, or evaporatively cooled in its presence. We investigate the dynamics 
using the quantum Boltzmann equation (QBE) for the Bose-Fermi gas mixture. The QBE has been applied to 
describe the population dynamics for a Bose gas by a number of authors who solved it by direct integration |p^| , |l4| , 
trajectory simulation [jT5|] , and Bird's simulation method ]T^| . The first two sets of authors make the assumption 
that the distribution is ergodic at all times |17| to simplify the numerical computation, whereas in Bird's method 
collision dynamics are constructed from simulated particle trajectories. In this paper we solve the QBE in the ergodic 
approximation by direct integration of the coupled differential equations which describe the Fermi and Bose gas 
distribution functions. 

We consider an harmonically trapped gas of Nf fermions with Fermi temperature given by Ep = ksTp = 
hu^Nf) 1 / 3 , which is in thermal contact with a similarly confined gas of 2V& bosons with condensation tempera- 
ture ksTc = huj(Nf ) /1.202) 1 / 3 . By evaporatively cooling the Bose gas through Tc down to a temperature rjTc 
(0 < r) < 1), the Fermi gas will equilibrate to the Fermi temperature provided Nf < if Nc /7 '.212, where Nc < Nf, is 
the number of condensed atoms remaining. Solution of the QBE gives both the energy and spatial distributions of 
the Bose and Fermi gases. An important limitation of our present treatment is the neglect of mean field effects on the 
dynamics. These are potentially important when the Bose and Fermi gas are degenerate and the mean field energy 
exceeds the trap frequency. In a recent paper Maimer [jl2| has used a simple mean field model to study the spatial 
distributions of a Bose-Fermi gas mixture at T = K . The distributions depend strongly on the relative sign and 
magnitude of the boson-boson and boson-fermion scattering lengths. A complete and numerically tractable approach 
to quantum kinetics which incorporates such mean field effects in addition to the population dynamics, is still under 



investigation. A qualitative discussion of the influence of mean field effects on our results is given in section IV. 

The remainder of this paper is organized as follows. In section |n| we discuss the equilibrium temperatures for an 
initially hot Fermi gas placed in thermal contact with a cold Bose gas. We also derive a simple dynamical model for the 
thcrmalization using the assumption that both gases are described by a Maxwell Boltzmann distribution at all times. 
In section III we describe a theoretical model based on the QBE for the Fermi-Bose mixture and the treatment of 



forced evaporative cooling. In section [Fy| we present results of solutions of the QBE which provide information on the 
kinetic temperature and spatial distributions of the mixture. Using recent atomic data we also present some results 
for the sympathetic cooling of the 40 K/ 39 K Fermi/Bose potassium isotopes. Finally, in section we summarize our 
conclusions. 



II. EQUILIBRIUM PROPERTIES AND COOLING WITHOUT FORCED EVAPORATION 

A. Equilibrium temperature 

If two gases at different initial temperatures are brought into contact they will rethermalize to a common equilibrium 
temperature which can be obtained using conservation of particle number and energy (over bar indicates average 
energy) 

etot =e/(0)+e 6 (0) 

= e>(oo) +e 6 (oo) (1) 

assuming that there are no losses during relaxation. For a Bose-Fermi gas mixture the equilibrium temperature 
can be found numerically using the Bosc-Einstcin and Fcrmi-Dirac distribution functions as follows. Given the number 
of bosons Nb, and fermions Nf, and the initial temperatures 21,(0), 2/(0) we set Tf = T b at each iteration step and 
compute the fugacity Zf for the Fermi gas by solving 

N f = J29(ei)(zJ 1 e^ kBT '+l)- 1 (2) 

with g{ti) the degeneracy of states with energy e,. We then compute the mean energy et(zf,Tf) for the Fermi gas 
and determine T b and z b for the Bose gas by the following iteration: set e b = etot — £f( z fjTf) and compute T b from 

e - 6 = ^e^)(C 1 e e4 / fc ^-l)- 1 . (3) 

Then compute by solving 

^^(e^-W^-l)- 1 (4) 

and repeat this iteration until [N b (z b ,T b ) - N b ]/N b < 1CT 7 and [e b (z b ,T b ) - e b ]/e b < 10~ 7 . We then set Tf = T b and 
repeat the procedure until (Tf - T b )/T f < 1CT 4 . 

For temperatures Tf > 0.5Tp and < Tq an approximate equation can be obtained using a classical Maxwell 
Boltzmann distribution for the Fermi gas and a Bose distribution with z b = 1 for the Bose gas. This yields 

„ m t T 4 (0) T/(0) 

with the Bose-Einstein function 54(1) ~ 1.082. If > Tc we can approximate the Bose-Einstein distribution with 
a Maxwell-Boltzmann distribution to obtain an explicit expression for the temperature 

_ NfTfjo) + mm 

lc °~ N f + N b ■ {b) 

In Fig. |l| we show the equilibrium temperature as a function of the initial Fermi gas temperature scaled in units of 
Tp. The Bose gas consists of 10 6 atoms at initial temperature T b (0) = O.lTp(O), where Tp depends on the number of 
fermions which is varied between 10 4 and 10 6 . In terms of the BEC temperature the initial temperature of the Bose 
gas lies in the range 0.02Tc < T b (0) < 0.2Tc- The full numerical results agree well with approximation Eq. (||) for 
Nf = 10 3 , 10 4 and 10 5 in the temperature range considered. This is shown by the curves at the bottom of the figure 
which almost overlap the approximate solution. The equilibrium temperature Too stays below the critical temperature 
Tc and in agreement with Eq. (||) T^/Tp does not depend on the number of fermions. If we increase the number of 
fermions or alternatively the initial Fermi gas temperature, the Bose gas will be heated above Tc and the equilibrium 
temperature shows the linear dependence on 2/(0) as derived in (^|). From our numerical studies we find that Fermi 
statistics become significant only for Tf < 0.5Tp (see also flgf] ), which in general can be reached with additional 
evaporative cooling of the gas mixture. 



B. Thermalization of a two-component mixture 



We consider Bose and Fermi gases brought into thermal contact in a confining potential, at different initial temper- 
atures T b (0) and 1/(0), respectively. We are primarily interested here in a regime prior to a stage of forced evaporative 
cooling, in which the Bose gas as well as the Fermi gas is non-degenerate. The thermalization may result in a signifi- 
cant alteration in the energy and spatial distribution of fcrmionic atoms in the trap. A simple dynamical description 
of the thermalization can be derived using the classical Boltzmann equation in the ergodic approximation. The mean 
energy of component k = (6, /) is given by 

e k (t) = J eF k (e,t) p(e)de, (7) 

where J- k (e,t) is the ergodic distribution function for component k, and p(e) = e 2 /2(Tno) 3 is the density of states for 
an harmonic trap of mean frequency to = (w I Wja) z ) 1/ ' 3 . For simplicity we asssume that both components see the same 
trapping potential. 

The time dependence of the energy of the Fermi gas is then obtained by integrating the Boltzmann equation |T^| 
over all energies 



dE f 1 1 
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dt To (hoj) 5 

[f f {E 4 ,t)F b {E 3 ,t) - Ff{E 2 ,t)F b {E u t)} , (8) 

where all energies written with upper case E are dimensionless, i.e, E = e/?uu, E m i n = min{_Ei, E 2 , E 3 , -E4}, and the 
natural timescale To is given in terms of the Bose-Fermi s-wave collision cross section a b f by 1/tq = (Twj) 2 ma b f / tt 2 Tl 3 . 
We have dropped the term which represents collisions between fermionic atoms assuming that the p-wave contribution 
is neglibible at the low energies under consideration. 

To get some qualitative insight into the thermalization dynamics we need to further simplify the model. We assume 
that the component distribution functions are Boltzmann-like at all times, and parametrized by time dependent 
fugacity Zk(t) and dimensionless temperature T k (t) = kBTk(t)/huj as follows, 

T[E r z k {t),T k {t)]=z k {t)e~ E l T ^\ (9) 
The average energy and particle number for a trapped Maxwell-Boltzmann distribution are given by 

E k = 7 ±n J dxx 3 e~\ (10) 

Nk = ^f 3 J dxx 2 e- x (11) 

from which it follows that z k — N k /T 3 , and E k = 3N k T k . The fugacity of component k can therefore be eliminated 
in terms of the particle number and mean energy. Another simplification which results from the approximation is 
that the scattering of two bosons does not alter the average energy of the Bose gas. 

It is convenient to write the average energy in terms of a dimensionless temperature, thus 

= / dx ± J dx z J dx z J dx 4 ^-x l S(x 1 +x 2 -x 3 - Xi ) {e~ Xi e- rXi - e^e-™ 1 ) 

= ^ P (T f /T b ), (12) 

§ WW, (13) 

where r = Tf/T bl and the rational function 

1 + 3r + 2r 2 - 2r 3 - 3r 4 - r 5 

P(r) = ; ■ (14) 

W (1 + r) 5 V ' 

The equilibrium temperature is obtained from dTf/dT b = —N b /Nf. Integrating from t = to 00, we recover 
equation (^|). The analysis shows that the timescale for relaxation is approximately T b To/Nf. In Fig. we compare 



the predictions of this simple dynamical model, with the full dynamics of the QBE as discussed in the following 
section. The solutions show that the time scale that describes equilibration is on the order of r r = T^To/Nf rj 0.043to, 
as predicted by Eqs. ( |l2| ) and (13). 

If forced evaporative cooling is applied to the gas mixture, atoms from the hot energy tail are removed. The rate 
of evaporation must be chosen to be much smaller than the relaxation rate of the gas as discussed above. Another 
way to determine the relaxation time is to calculate the initial mean collision rate 7 using the Boltzmann collision 
integral. Consider the collisions of "test" atoms, with energies E\, with atoms of energy, E2, into final states with 
energies E3 and E4. One sums over all initial energies E\, and divides by the number of particles to get the energy 
averaged rate per particle 
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in qualitative agreement with the simple model above. We note, however, that the relaxation rate we have discussed 
here should not be regarded as the characteristic time scale for condensation [p~9|—^T]| . 

III. QUANTUM BOLTZMANN EQUATION AND FORCED EVAPORATION 

A. Quantum Boltzmann equation 

The QBE for a harmonically confined Bose gas has been discussed elsewhere within the ergodic approximation 
|l5| , p2|| . The QBE for an interacting two component Bose- Fermi mixture in the harmonic trap can be written (t = t/ro) 

9{Ei)—^- — a b ^2 SE i +E j ,E h ,B t g(Ei,Ej,Ek,Ei) 

Ej ,E k ,Ei 

[b Ek b El (l + b Ej )(l + b Ei ) - b Et b Ej (l + b Ek )(l + b El )] + 

^ ^Ei+Ej,E h +Eig(Ei, Ej, Efc, Ei) 
Ej,E k ,Ei 

[b Ek f Ei a - f E] )(l + b El )-b E J E] (l - f El )(l + b El )] , (16) 

df 

9i R i)-^-^ a f S E l +E,,E k +E l g(Ei,Ej,E k ,Ei) 

Ej,E k ,Ei 

[fE k f El (l ~ f Ej )(l - f Ei ) - f Ei f Ej (l - f Ek )(l - f El )] + 

Y S E t +E j ,E h +E l g(Ei,Ej,E k ,Ei) 

Ej,E k ,E t 

[b E J Ek (l - f Ei )(l + b E ,)-b Ej f Ei (l - f Ek )(l + b El )} , (17) 

where b n and /„ are the number of Bose and Fermi atoms in state E n . The collision matrix elements are approximately 
given by 

g{E i ,E j ,E k ,E l )=g{E mixi ), (18) 

with g(E n ) the degeneracy of energy level E n , and -Emin is the minimum energy of all four energies involved in the 
scattering process, as defined earlier. Although this approximation is not quantitatively accurate for the lowest several 
states of the trap, it is sufficient to illustrate the main qualitative features of sympathetic cooling. In an isotropic trap 
the degeneracy of the energy state e n = Tiw(n — 1), n = 1, 2 • • •, is g{E n ) — n(n + l)/2. The coefficients a b — (Jbb/^bf 
and Q.f = Gff/crbf, give the ratios of the cross sections for boson-boson and fermion-fermion scattering, respectively, 
to the boson- fermion cross section. Exchange symmetry leads to af — since the s-wave cross section vanishes for 
identical fcrmions. Of course, this is the reason we must employ sympathetic cooling with a Bose gas refrigerant. 



B. Forced evaporative cooling 



Evaporative cooling in magnetic traps is performed by inducing transitions to untrapped states with a radio- 
frequency field. This is modeled here by the following procedure. Particles that are scattered into states with energy 
larger than the time varying cut-off energy E cu t(T) are lost. The latter is a given decreasing function of time in 
the case of forced evaporative cooling. A particle may be scattered into a state with energy Ei < E cut (T), by two- 
body scattering of atoms in states with energies Ek and Ei, i.e., Ek, Ei — ► JSj, Ej, in which Ej > E cut (T), so that 
one particle is lost from the trap. Similarly a particle in energy level Ei can be scattered out of this energy level, 
Ei, Ej — > Ek, Ei, resulting in one particle loss from the trap when either Ek > E cu t(r), or Ei > _E cut (r). Explicitly, 

(a) The gain process for energy level Ei 

dn Ej<2E out (T) 

g{Ei)—^ L = a ^ SE^+E^E^EigiE^Ej^kiE^nE^Eiil ±n B J- (19) 

E j >E cut (r)>E k ,E l 

(b) The loss process for energy level Ei 

dn E fc <2B cut (r) 

9( E i)—r J -=-2® X! 5 Ei +E J .E k .E l g{E l ,E j ,Ek,Ei)nE t nE J {l±n El ), (20) 

E k >E mt (T)>E,,Ei 

where nE t denotes the distribution function for Fermi or Bose atoms with energy Ei, as appropriate. 

The kinetics of forced evaporative cooling is modeled by adding these terms to the QBE, Eqs. (|l^) and (|l^), which 
include all two-body collision processes between initial and final states i,j,k, and / with energies Ei, Ej, Ek, Ei < 
E cu t(r) conserving the total number of particles in the trap. 



IV. RESULTS AND DISCUSSION 



In this section we illustrate the dynamics of sympathetic cooling of Bose-Fermi gas mixtures, through their energy, 
state and spatial distribution functions. The spatial distributions of confined degenerate Bose and Fermi atomic gases 
are q uite diffe rent. An ideal Bose condensate has a size determined by the quantum width of the trap ground state 
I = \Jh/2Muj, whereas the size of a Fermi gas is governed by the Fermi width Rp = (Ep/2Mu; 2 ) 1 ^ 2 , which scales as 

1 /6 

Rp ~ NJ I, as a result of the Pauli exclusion principle. For an interacting Bose gas, with positive scattering length, 
the condensate is larger than I, and for strongly condensed gases its size can be estimated using mean field theory 
in the Thomas- Fermi approximation |2^| . Mean field effects, which can be significant well below the condensation 
temperature, are not included in our model. These may be important in the final stages of cooling if the Bose gas 
is already strongly condensed at this stage. Our illustrations of the spatial distributions of both Fermi and Bose gas 
employ the universal scaling described by Butts and Rokhsar , who showed that for a harmonically trapped ideal 
Fermi gas at T = 



/ > N f 8 

nf{r) = W F V 2 



\Rf) 



3/2 



(21) 



It should be remembered that with evaporative cooling the number of fermionic and bosonic atoms is a time dependent 
variable, and therefore so is Rp. In the figures we always scale with respect to the instantaneous value of Rp. 

In Figs. ^ and || we present the rethermalization of a non-degenerate Fermi gas immediately after it is placed 
in thermal contact with a Bose gas which is initially at the Bose condensation temperature Tq- The calculation is 
performed by numerical integration of the QBE without any forced evaporative cooling. In Fig. || the temperature 
of the Fermi gas alters considerably over a timescale ToTf/Nb [Eq. ([l2|)]. Initially a hot tail of atoms extends to the 
trap extremities and during the early stages of equilibration the Fermi gas distribution deviates significantly from 
a Fermi-Dirac distribution which is fit to the average energy and particle number. The gas then equilibrates to a 
non-degenerate state as can be seen from inspection of the peak of the spatial distribution function, n/(r = 0) <C 1, 
|l8|| . The Bose gas has one hundred times more particles than the Fermi gas, and completely envelopes the Fermi gas 
at all times. Fig. || compares the simple model of thermalization discussed in the last section with the QBE. The 



model is very good in the early stages, but the agreement deteriorates in the intermediate regime before steady state 
is achieved. 

In Figs. ^ and we consider the forced evaporative cooling of both gases. In contrast to Figs. ^ and |^ there are 
10 5 bosons and 10* fermions, initially. The forced cooling begins after the initial equilibration stage during which the 
boson temperature increases (Fig. |5|), following thermal contact of the two gases. The evaporative cooling time scale 
is chosen to be To, which is much longer than both the relaxation time scale of the one component Bose gas [Eq. jl^)] 
and the relaxation time scale of the Fermi gas with the Bose gas [Eq. (f[2|)] . The Bose gas energy distribution shows 
the formation of the condensate and the corresponding spatial distribution contracts to that of the condensate with 
a small thermal component. The degenerate Fermi gas is then exposed and its spatial distribution is close to the 
zero temperature limit [Eq. ([2l|)] which has a maximum density at the trap center rif(0)R F /Nf — 8/ir 2 « 0.81. The 
inset shows the state occupancy for the Fermi distribution with the characteristic smearing of the Fermi surface at 
finite temperature, and near unit occupancy for low lying levels. It is interesting to note that evaporative cooling of 
the Fermi gas still proceeds at later times when the spatial overlap between fermions and bosons is mainly in a small 
region at the center of the trap where the condensate is located. The collisions which result in cooling involve orbits 
of hot fermions through the trap center where they collide with cold condensed bosons. At this stage evaporation 
mainly results in depletion of fermions as can be seen in Fig. |^. We also simulate the case when the evaporative 
cooling involves only the loss of Bose atoms from the trap. In Figs. || and ^ we consider the same initial conditions 
as for Figs. [| and |5|, but only ramp down the cut-off energy for the bosons. The results are qualitatively the same 
as in the former case where both Bose and Fermi gas particles evaporate, except that we end up with more particles 
left in the Fermi gas. As mentioned earlier we have not included the effect of the bosonic mean field which can alter 
the spatial distributions of each component depending on the mean field strength, the ratio of the scattering lengths, 
and the particle numbers of both components fl2]| . 

A possible experimental scenario for sympathetic cooling of a Bose-Fermi mixture involves two isotopes of potassium. 
Recent calculations predict the s-wave scattering length for the bosonic isotope 39 K to be a( 39 K) = a& = 4.3 (nm) 
with corresponding cross section abb = 87r 2 a 2 and the s-wave scattering length between 39 K and the fermionic isotope 
40 K to be a( 40 K - 39 K) = a bf = 2.5 (nm) with cross section a bf = 4 7 r 2 ag / @. Using A/( 40 K) = 6.6 x IQ-^kg and 
u! = 400 (Hz) sets the time scale To = 1254 (s) « 20(min). Mean field effects become important when the mean field 
strength -E moa n = 47r7i 2 afc/(Af( 40 K)fio) is of the order of the level spacing htu 25|. At the onset of BEC the density 



profile of an ideal Bose gas is almost Gaussian and the mean density of the ground state becomes no = Nq/tt 3 ^ 2 1 3 
The ratio of the mean field strength to the trap energy is then 

^mean / 1 &6 AT - ^ n . . n — 3 1 



= = , = x 10 -d JVo< (22) 

TlUJ V 7T I 

Further discussions of the influence of mean-field effects on our results is given below. In Figs. J| and ^ we present 
an example of evaporative cooling for a spin polarized K— 39 K mixture of 10 6 bosons and 10 fermions at initial 
temperatures T b = T c = 0.3 (^Kjand T f = 7.2T F = 1.8 (/xK) and cut-off energy E cut (0) = 1000. We estimate the 
boson scattering rate using Eq. ( |l5| ) which yields w 5 x 10 3 /to. From r = until r = 0.04 the cut-off energy remains 
at E cu t = 1000. During the thermalization the Bose distribution completely envelopes the Fermi distribution which 
deviates significantly from its equilibrium distribution. As one can see from Fig. || the number of fermions almost 
remains the same whereas a large number of bosons are evaporated. After r = 0.04 the cut-off energy is ramped down 
exponentially with rate 7 OV ap = 100/to until r = 0.064 = 80.25 (s). At the early stages of the forced evaporation the 
number of bosons decreases until most of the bosons are in the condensate and evaporation mainly leads to depletion 
of fermions. The simulation shows that the Fermi gas can be cooled to a temperature Tf ~ O.lTp with more than 
2 x 10 4 fermions left in the trap. 

In practice there are some additional issues that must be considered. The isotopes have different mass and magnetic 
moment, this means that in general the clouds will be displaced with respect to one another due to the combined 
effects of gravity and the magnetic trapping force || . Sympathetic cooling can only proceed efficiently if good overlap 
between the gases is maintained p6[ . Even if we assume this has been achieved, the difference in magnetic moments 
will cause the trap frequencies to be different for the Bose and Fermi gas. For example if the bosonic 39 K ( / = 7/2, 
where I is the nuclear spin) is polarized in the state \F = 2, Mp = 2), and the fermion isotope 40 K (I = 4) is polarized 
in the state \F — 7/2, Mp — —7/2), the trap frequencies would be in the ratio of 7 : 9. In our calculation we assume 
the trap frequencies and masses to be identical. 

For degenerate Bose and Fermi gases mean field effects can strongly influence the spatial distributions. Here we 
discuss how these effects qualitatively change the stationary distributions presented above, using the zero temperature 
model discussed by M0lmer p2| . If the number of bosons is much larger than the number of fermions then as abf 
increases relative to Of, (abf, a& > 0), the fermion distribution is displaced further and further outside of the central 
core of the trap occupied by the bosons. When the particle numbers are similar the bosons are displaced outside the 



fermion region in the same limit. On the other hand, for large negative fermion-boson scattering length the stationary 
distributions may be unstable and dynamically change as a result. In the latter case mean field effects qualitatively 
change the nature of the problem, and thus our results will not apply. 

In our simulations the number of bosons is always much larger than the number of fermions. In this case we can 
to first order neglect the influence of the fermions on the the boson spatial distribution. The boson density is then 
given by the Thomas-Fermi approximation 



n b (r) = [fi - V ext (r)]/a b 



(23) 



where I4 x t(r) = Muj 2 r 2 /2 denotes the harmonic trapping potential, and the chemical potential /i is fixed through the 
condition J drnb{r) — N b . Explicitly, this yields 



M=[( 



raw 15 



2/5 



(24) 



The bosonic mean-field produces the well known broadening of the boson density distribution relative to the ground 
state length of the trap. The corresponding spatial distribution of fermions can be found from the equation fl2|] 



2M 



[6^ 2 n/(f)] 2 / 3 + (1 - ^-)F ext (f) + ^L^ = E F 



where Ef is determined by J dfn / (r ) ~ Nf . Explicitly the density is given by 



Nf 8 



1 _ 06/ jj_ _ _ 06/^ / _r_ 
en, Ep a,b \Rf 



-i 3/2 



(25) 



(26) 



which may be compared with equation ( |2l|) for an ideal gas. Clearly the mean field effects cause a broadening in 
the spatial distribution of the fermion cloud for a&/ < a^. If a&/ > a,b the fermions experience an inverted harmonic 
oscillator potential near the origin which repells them from this region. In our results presented above, and in 
particular for the example of the potassium isotopes, we always have abf < Q-b/%- As a result mean field effects will 
result in a broadening of both Bosc and Fermi gas spatial distributions, but not a relative displacement of the clouds. 



V. CONCLUSION 



We have discussed the cooling of a confined non-degenerate Fermi gas to quantum degeneracy using an ultracold 
Bose gas coolant and evaporative cooling. Results for the stationary distributions and dynamics based on solutions of 
coupled QBE equations for the Bose- Fermi mixture were presented. These include investigations of the use of forced 
evaporative cooling to enhance the degeneracy of the Fermi gas. While the QBE does not include mean field effects, 
which are potentially important in the quantum degenerate regime, we have discussed their qualitative effects on the 
results presented, when a& > a&/ > 0. In this instance mean fields lead to a broadening of both the Bosc and Fermi 
spatial distributions, but not a relative displacement of the clouds. 
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FIG. 1. Equilibrium temperature as a function of the initial Fermi gas temperature in units of Tp. The number of fermions 
is Nf = 10 3 (-), Nf = 10 4 (:) JV> = 10 5 (.), N f = 2 x 10 5 (o), Nf = 4 x 10 5 (0) , iV> = 10 6 (+). The solid lines show the 
results using Eq. (H) and Eq. (H) respectively. The number of bosons is iV;, = 10 6 and the initial temperature is always chosen 
as 51,(0) = 0.15V corresponding to f b (0) = 0.0197c, 0.042Tb, 0.09Tb, 0.135b, 0.142Tb, and 0.19Tb- 

FIG. 2. The temperature of the Bose gas (solid line) and of the Fermi gas (broken line) as a function of time. The dotted 
line denotes the temperature obtained by solving the differential equations assuming a Maxwell Boltzmann distributions with 
time dependend fugacity and temperature. 

FIG. 3. The Fermi (solid line) and Bose (dotted line) distribution functions at times r = 0, 0.003 and 0.016. The graphs on 
the left show the number of atoms as a function of energy, the graphs on the right the spatial distribution. Initial conditions 
are N b (0) = 10 5 , N f (0) = 10 3 , and 51,(0) =f c = 43.7, 5>(0) = 55V = 81.8. Both gases are in contact at r = and then relax 
to the equilibrium temperature Too = 44. The broken lines denote the fit to the Fermi distribution and coincide with full lines 
at tau = and r — 0.016. 

FIG. 4. The Fermi (solid line) and Bose (dotted line) distribution functions at times r = 0, 0.9, and 2.0. The graphs 
on the left show the number of atoms as a function of energy, the graphs on the right the spatial distribution, the inset the 
number of fermions per energy divided by the degeneracy. Initial conditions are Nb(0) = 10 5 , Nf(0) = 10 4 , 51,(0) = 5b = 43.7, 
5>(0) = 55V = 186. From r = until r = 0.04 the cut-off for the Bose gas is E cut = 500 and for the Fermi gas E cut = 1000. 
After r = 0.04 the cut-off energy is ramped down with a rate 7 = 1.0 for both gases starting at Eo = 500. [E cu t(r) — e~ 1T Eo]. 
The fit to the Fermi distribution is also drawn as a broken line. The temperatures from top to bottom are Tf(r) = 55V, 0.65V 
, 0.145V, and T b (r) = T c , 0.5Tb, 0.1Tb 



FIG. 5. Upper graph: The number of bosons and fermions remaining in the trap. Lower graph: The temperature of the 
gases as a function of time. The dotted lines show the critical temperature for the Bose gas and the Fermi temperature. 



FIG. 6. The Fermi (solid line) and Bose (dotted line) distribution functions at times r = 0, 0.9 and 2.2. The graphs on the 
left show the number of atoms as a function of energy, the graphs on the right the spatial distribution, the inset the number 
of fermions per energy divided by the degeneracy. Initial conditions are Nt(0) = 10 5 , Nf(0) = 10 4 , 51(0) = To = 43.7, 
5/(0) = 55V = 186. From r = until r = 0.04 the cut-off for the Bose gas is E bcut = 500 and for the Fermi gas 
Ef cu t — 1000. After r = 0.04 the cut-off energy is ramped down with a rate 7 = 1.0 only for the Bose gas starting at 
i?o = 500. [E cu t(r) = e~ JT Eq]. The fit to the Fermi distribution is also drawn as a broken line. The temperatures from top to 
bottom are 5>(r) = 55V, 0.55V, 0.145V, and T b (r) = 5b, 0.5Tb, 0.145b- 

FIG. 7. Upper graph: The number of bosons and fermions remaining in the trap. Lower graph: The temperature of the 
gases as a function of time. The dashed lines show the critical temperature for the Bose gas and the Fermi temperature. 

FIG. 8. Sympathetic cooling of the fermionic Potassium isotope 40 K by the bosonic isotope 39 K. The Fermi (solid line) and 
Bose (dotted line) distribution functions at times r = 0, 3.6 x 10~ 4 , 4.7 x 10 -2 , and 6.1 x 10~ 2 . The graphs on the left show the 
number of atoms as a function of energy, the graphs on the right the spatial distribution, the inset the number of fermions per 
energy divided by the degeneracy. Initial conditions are N b (0) = 10 6 , N f (0) = 10 5 , 51(0) = Tb = 94.1, 5>(0) = 7.2T> = 590.4. 
From t — until r = 0.04 the cut-off for both gases E cut = 1000. After r = 0.04 the cut-off energy is ramped down with a 
rate 7 = 100.0. The broken line denotes the fit to the Fermi distribution . The Fermi gas temperatures from top to bottom 
are 5>(t) = 7.25V, 1.625V, 0.675V, and 0.145V. The Bose temperature is f b (r) = 5b, 1.02Tb, 0.63Tb, and 0.125b- 

FIG. 9. Upper graph: The number of bosons and fermions remaining in the trap. Lower graph: The temperature of the 
gases as a function of time. The dashed lines show the critical temperature for the Bose gas and the Fermi temperature. 
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